Absence of Thermalization in Nonintegrable Systems 
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The question of how quantum many-body systems in 
nonequilibrium eventually equilibrate and assume properties 
resembling the ones familiar from statistical mechanics has — 
quite unsurprisingly — a very long tradition (T). In closed 
systems not all observables can equilibrate. However, it is 
generally expected that in sufficiently complicated quantum 
many-body systems at least some physically relevant quanti- 
ties should seemingly relax to equilibrium values. Recently 
this old question has received an enormous amount of atten- 
tion and there have been significant new insights B ffl"8l . 

This renewed attention is partly driven by new mathemat- 
ical methods becoming available, partly by novel numerical 
techniques, and in parts by experiments that make it pos- 
sible to probe coherent nonequilibrium dynamics under the 
controlled conditions offered by cold atoms in optical lat- 
tices lfl9l . Theoretically, among other approaches, the ques- 
tion of how quantum many-body systems relax locally has 
been investigated in the light of the "eigenstate thermaliza- 
tion hypothesis" (ETH) ||2] |3J, quantum central limit theo- 
rems 0, Anderson localization [7|, dynamical instances of 
concentration of measure arguments or ideas of relaxation via 
dephasing [8-13, 16 1, and numerically using time-dependent 
density-matrix renormalization group (DMRG) |[2D . Despite 
this enormous effort, major questions remain open and the ex- 
isting results do not yet draw a coherent picture. What seems 
to have become consensus, however, is that the following ex- 
pectation holds true: Nonintegrable systems thermalize. 

In this letter we show that generally, this is not quite true. 
We do so by establishing a link between the entanglement 
in the eigenbasis of a quantum many-body system with what 
could be called the thermalization potential of the system. We 
will investigate situations in which systems equilibrate, in the 
sense that all local observables will be close to some equilib- 
rium value at most times, but those values turn out to depend 
on the details of the initial state. This general rigorous state- 
ment is exemplified numerically by studying a small natural 
nonintegrable XYZ-type spin chain model. In previous ap- 
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Figure 1 . (color online) A many-body quantum system consisting of 
a small subsystem S, and the remainder B, which acts as a "bath." 



proaches (e.g., in Ref. Il23l ). similar complementing observa- 
tions have been made by simulating the model's time evolu- 
tion explicitly. However, such simulations can only trace the 
system's behavior for a finite amount of time and become un- 
reliable for long times. Our analytic results have applications 
far beyond this particular model: they yield general condi- 
tions for the absence of thermalization. This gives a natural 
counterpart of the ETH, and it relates the thermalization of 
isolated quantum systems to the presence of entanglement in 
the energy eigenbasis. 

Setup and notation. When using terms and concepts bor- 
rowed from classical statistical mechanics, such as ergodicity, 
equilibration, thermalization, initial state independence, and 
integrability, we aim at being careful and precise whenever 
referring to one of these terms. We work in the pure state 
quantum statistical mechanics model with a system and bath 
setup with a global pure state and unitary time evolution flUE]- 
[121 . We are mostly interested in the case where the full system 
is composed of many interacting small systems and the sub- 
system corresponds to a small subset of sites and the bath is 
simply the remainder.In such systems the individual subsys- 
tems act as baths for each other and the collective dynamics 
can lead to self-thermalization of the whole system. To be spe- 
cific, we will consider arbitrary quantum systems equipped 
with a Hilbert space H of finite dimension d that can be di- 
vided into at least two parts, i.e., H = Hs <S) Hb, which we 
will call the subsystem S and the bath B, and which are de- 
scribed by Hilbert spaces of dimensions ds/B = dim(Ws/s). 
We assume that at every time t the joint system is in a pure 



2 



state ipt — \if>t){ l f , t\> evolving unitarily. The reduced states on 
subsystem and bath are denoted using superscript letters, such 
as ip s = Tr# [?/>]. We denote the Hamiltonian of the full sys- 
tem by Jtf and its eigenvectors and eigenvalues by \Ef.) and 
Ek, k = 1, . . . , d (regardless of degeneracies). 

Thermalization. Thermalization is a complicated process. 
For it to happen a system must exhibit certain properties, each 
of which captures a specific aspect of thermalization. It is 
instructive to consider each of them separately ll9l [TD : 

1) Equilibration: The tendency to evolve towards equilib- 
rium is a key assumption in classical statistical physics and 
part of the second law of thermodynamics. In contrast to that, 
in the framework of pure state quantum statistical mechanics, 
equilibration for almost all times can instead be proven to be 
a consequence of unitary time evolution |9|. 

The remaining conditions specify what properties the equi- 
librium state of a small subsystem should have and of which 
of the initial conditions it should be independent. 

2) Subsystem initial state independence: The equilibrium 
state of a small subsystem should be independent of the initial 
state of that subsystem. This aspect of thermalization will be 
the main subject of the present work. 

3) Diagonal form of the subsystem equilibrium state: The 
equilibrium state of a small subsystem should be (close to) 
diagonal in the energy eigenbasis of its self-Hamiltonian (8). 

4) Bath state independence: It is generally expected that 
the equilibrium expectation values of local observables on a 
small subsystem are almost independent of the details of the 
initial state of the rest of the system, but rather only depend 
on its macroscopic properties, such as the energy density. 

5) Gibbs state: Ultimately, one would like to recover the 
standard assumption of classical statistical physics that the 
equilibrium state is of (or at least close to) a Gibbs state 
lj s Rj e - ^ Jifs with an inverse temperature (3 and a self 
Hamiltonian Jtffs- As in conventional statistical physics this 
can only be expected to be true if the coupling is weak but 
nonperturbative and the bath has a spectrum that gets expo- 
nentially dense for higher energies lfT4l[T5l . 

Integrability. In classical mechanics integrability is a 
well-defined concept [27 1. In quantum mechanics, despite the 
common use of the term "integrable," the situation is much 
less clear 11221 . and different criteria are being applied in the 
literature. The most common notions of integrability are the 
following: 

(A) There exist n independent (local) l28l conserved mutu- 
ally commuting linearly independent operators, where n is the 
number of degrees of freedom (see, e.g., Ref. 11241 ). In con- 
trast to the classical situation 1271 . this does not necessarily 
imply that the system is "exactly solvable." 

(B) Identical with criterion (A), but with linear indepen- 
dence replaced by algebraic independence E4l . 

(C) The system is integrable by the Bethe ansatz [24|. 

(D) The system exhibits nondiffractive scattering [24|. 

(E) The quantum many-body system is exactly solvable in 
any way. Of course, this criterion is subject to the ambiguity 
of a lack of imagination of solving a given model. 



Equilibration. Quantum mechanics of closed systems is 
time reversal invariant and thus equilibration in the usual 
sense is impossible. Therefore we use an extended notion of 
equilibration and say that an observable A equilibrates if its 
expectation value Tr[vlV>t] is close to some value for almost 
all times t. Of particular interest are local observables, i.e., 
observables that are sums of terms that each act only on small 
subsystems. Saying that all observables on some subsystem S 
equilibrate is equivalent to saying that the state ipf of the sub- 
system equilibrates, by which we mean that there exists a state 
p° such that ipf is almost always physically indistinguishable 
from p s ', in the sense that their trace distance T>(ipf, p s ) — 
\\\ipf - p s \\i = max o <A<iTr[A'0f] -Ti[Ap s ] is small for 
almost all times t. If the expectation value Ty[A ipt] of an ob- 
servable A equilibrates in the sense defined above, then it must 
equilibrate towards its time average Tr[A-0 t ] = Tr[A?/> t ]. 
This is an obvious but important observation. A good under- 
standing of the properties of the time averaged state 
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ui = ipt = lim ~ I ipt at ( 1 ) 
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is thus key to understanding equilibrium properties. 

Generalized Gibbs ensemble. Every Hamiltonian Jf? de- 
fines a set of conserved observables. In the nondegener- 
ate case they are exactly the linear combinations of projec- 
tors onto the eigenstates of J^, in the degenerate case they 
are the observables with support on the blocks correspond- 
ing to the degenerate subspaces. Clearly, the time average 
u) = ipt of the state i/jq itself is given by ui = P(ipo), where 
P{^o) = J2j njipQir-j, where 7Tj = Ylkel, \ E k)(E k \ are the 
projections onto (possibly degenerate) eigenspaces, Ej~ = E\ 
for k,l E Ij. Every state p that gives the same values for 
all conserved observables as tpQ satisfies P(ipa) — P(p) and 
uj = P(tpo) is the state having maximum entropy among all 
such states. This follows directly from the pinching inequality 
(Theorem V.2.1 in Ref. ||20l ) since the von Neumann entropy 
is Schur concave. All equilibrium expectation values can be 
calculated from the maximum entropy state u. This is a quan- 
tum version ofJaynes '-principle and was recently conjectured 
as generalized Gibbs ensemble in Ref. 1 18 1. 

Moreover, under the assumption of nondegenerate energy 
gaps it can be rigorously proven under which conditions equi- 
libration (but not necessarily thermalization) happens ll9l [T3l . 
The certificate quantifying the quality of equilibration is the 
effective dimension of the time averaged state d cS (u) = 
1/Tr[w 2 ], which, for quenches to nondegenerate Hamilto- 
nians, is identical to the inverse of the time average of the 
Loschmidt echo and the inverse participation ratio (IPR) of 
the initial state |26|. The main result of Ref. |9| is 



and we call C eq (tpo) the equilibration coefficient of the initial 
state ipo as it bounds the trace norm equilibration radius. 
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Main result. In systems that behave thermodynamically 
the equilibrium expectation values of local observables on 
small subsystems should be independent of the initial state 
of the subsystem. A previous positive result in this direction 
was made in Ref. |9] (see also Refs. HUE]). Here we follow 
a converse approach and proof a sufficient condition for the 
absence of initial state independence. 

A quantity that will play an important role in our main re- 
sult is the effective entanglement in the eigenbasis, given for a 
nondegenerate by 
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c k \ 2 V(Tr B \E k )(E k \,il>(j] 



(3) 



with Cfe = (Ek\tpo)- This quantity is small, if most energy 
eigenstates either resemble locally the system's initial state 
■00 > or are globally almost orthogonal to i/jq. As will become 
apparent later, this is in particular the case if the reductions 
of the \E k )(E k \ are close to a basis for S. This can be in- 
terpreted as a natural counterpart of the ETH [12 j: If "most" 
energy eigenstates have reduced states close to some p s , then 
the system will relax locally to p . 

We will now show that a small value of R implies that initial 
state independence is not satisfied. Remarkably, this is not a 
matter of time scales: It will not only take a long time to relax; 
but one will rather encounter a memory for almost all times. 



Theorem 1 (Nonthermalization). The physical distinguisha- 
bility of the two local time averaged states uj s ^ andto s ^ of 
two pure initial product states ipQ = ® <i>a^ % \ * € 

{1,2} evolving under a nondegenerate Hamiltonian Jff is 
large in the sense that 



In the degenerate case, the quantity R has to be replaced by 

k 



That is to say, subsystems remain distinguishable if they are 
initially well distinguishable and one has little effective en- 
tanglement in the eigenbasis. Note also that the environment 



states 



and 



iB(2) 



can be taken to be identical. 



Proof. If the Hamiltonian ,J%° is nondegenerate, dj s ^ = 



\$\*Tr B \E k )(E k \ 



and thus 
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c k \ z Tv B \E k )(E k \ 
<$>k| a ~K -Tr B \E k )(E h \ 



The desired result then follows from T>(ip ,tp ) 



< 



S(l) 



-V(uj 



S(l) ,,,S{2)\ 



-V{w s{ - 2 \% (2) ). In the 



degenerate case, the same argument can be followed for the 
projectors n k onto the respective eigenspaces. □ 



The intuition that R will be small when the \E k ) have cer- 
tain properties can be made rigorous in the case where the 
Tr# \E k )(E k \ are close to a basis of S. In this case we can 
show that there exist many initial states for the bath that lead 
to a small R and a large effective dimension at the same time, 
thus causing "equilibration without thermalization." This is 
shown using Haar-measure averages, from which the exis- 
tence follows ll29l . 

Theorem 2 (Entanglement in eigenbasis). For every or- 
thonormal basis {\i)} for S and every initial product state 
with ipQ — for some i and with Haar random initial 

bath part <p B , the effective entanglement in the eigenbasis for 
nondegenerate Jrff is on average upper bounded by 

E^bR(^®^) <25d s 

where 5 = max fe 5 k with S k = min^ V(Tr B \E k )(E k \, 
being the geometric measure of entanglement of the eigenstate 
\E k ) with respect to the basis {\i)}. 

Proof. Note that Tr[Tr B \E k ) (E k \ipg] < 1 - 
T>(Tr B \E k ){E k \,ip$) 2 and thus all nonzero |c fc | 2 in Eq. 
<|3j can be upper bounded by 

Ti[\E k }(E k \(^^^ 



^l(l-V(Tr B \E k )(E k \,ip s ) 2 ). 



Tr[Tr B \E k )(E k \^ 

As (1 - V(Tr B \E k )(E k \,^) 2 )V(Tr B \E k )(E k \,^) < 
2 6, we have 



k 



Tr[\E k )(E k \(^®<f> B )] 
Ti[Tr B \E k )(E k \ipS] 



Averaging over all pure states <p B gives the mean 1 /d B and 
the sum in the last line is thus upper bounded by dg. □ 

Note that Theorem|2]implies that whenever a basis { \i) } ex- 
ists for which 6 is small, then for every i there exist many bath 
states 4> B such that R(\i)(i\ ® 4>q ) < 2Sd s (29). Further- 
more, almost all of them will lead to a high effective dimen- 
sion J9HT2). Obviously the argument can be further strength- 
ened by maximizing 6 only over a subspace that contains most 
of the probability weight of t/j n : This allows some of the 8 k to 
be large, as long as the corresponding \c k \ 2 are small. 

Application to a nonintegrable model. The model we con- 
sider is a spin- 1/2 XYZ chain with n sites with random cou- 
pling and on-site field. The Hamiltonian is 



i '//^ 



(4) 



where of N = (of af +l ,aj o^.of af +l ) T in terms of the 
Pauli matrices at site i, and hi and the components of hi are 
i.i.d. normal distributed random variables with zero mean and 
standard deviations oo and <j\, respectively. The model is 
closely related to the one studied in Ref. [7]. With unit prob- 
ability the Hamiltonian is nondegenerate and has nondegen- 
erate gaps. We investigate the equilibration properties of the 
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proximately conserved quantities can exist even in noninte- 
grable models. Such models may not saturate Lieb-Robinson 
bounds, i.e., there probably is no ballistic propagation of in- 
formation. Certainly, interesting physical candidates for such 
models are to be found in disordered systems: The Anderson 
model for example has eigenfunctions that are exponentially 
localized with high probability ||251 . It is the hope that this 
work stimulates further research on this connection. 
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Figure 2. (color online) The subsystem is taken to be the first site 
S = 1 in the spin chain of n sites, other choices give qualitatively 
the same results. For each of the product eigenvectors \E®} of Jifo 
we compare the equilibration properties of tp^ — \E%) with that of 
t/)q 2 ' = o§\E®} (i.e., the same state but with the first spin flipped) 
under the dynamics of Jif with <7i = 0.4 ao. Panels (a), (b) dis- 
play averages over eigenstates: (a) Average geometric measure of 
entanglement K(5k) with respect to the Jif o eigenbasis and average 
distance of the reduced time averag ed states E(£>(w S(1) , w s(2) )). 
(b) Average effective dimension and equilibration coefficient [see 
Eq. {2)]. Panels (c), (d) show quantities optimized over eigen- 
states: (c) Maximum distinguishability maxt A(\E^)) where 
A(|4°>» = 2>( W s «, W s < a >) _ c^W) _ Ccq (4 2) ) (A > 
ensures distinguishability for almost all times. See the inset for an 
artists impression.), (d) Effective dimension and equilibration coef- 
ficient of the state maximizing A(\E^}). All quantities have been 
averaged over 100 samples from the random Hamiltonian ensemble 
Q. The error bars represent the standard deviation. d e increases 
rapidly with n: hence, equilibration gets better, while the time av- 
eraged states on S remain well distinguishable. Remarkably the ob- 
servable that best distinguishes the equilibrium states is a§ . We find 
"equilibration without fhermalization" for a very natural observable. 

eigenstates of Jj? after a quench to via exact diagonal- 
ization for <j\ — 0.4 cto. Hence the integrability breaking 
term is far from being a small perturbation and ^ is 
nonintegrable in the sense of all of the aforementioned defi- 
nitions of integrability. According to the widely accepted be- 
lief E|[2l][23|, one would therefore expect to find thermaliza- 
tion. However, the numerics suggests that after the quench all 
local observables equilibrate, but retain memory on the ini- 
tial conditions and thus initial state independence is violated. 
This conclusion is reached not by keeping track of time evolu- 
tion, but rather by checking the conditions of Theorem[T](see 
Fig. |2j. It is a challenge to construct nonintegrable models 
without disorder that violate initial state independence. 

Summary and conclusions. We have established rigorous 
results that identified a lack of entanglement in the energy 
eigenbasis as the reason for an "equilibration without ther- 
malization" phenomenon: all local observables equilibrate but 
retain memory on their initial values for infinitely long time. 
By considering a particular model we exemplify that such ap- 
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EPAPS 

Connection between effective dimension, Loschmidt echo and 
inverse participation ratio 

The effective dimension of the time averaged state, the 
time average of the Loschmidt echo and the inverse partici- 
pation ratio of the initial state are closely related in the non- 
degenerate case. The inverse participation ratio (IPR) of a 
state V>o 1S defined to be 



IPR(Vo) = 



Ck\ 



where Ck = (Ek\i^o) are the overlaps with the energy eigen- 
vectors {|£?fc)} of the Hamiltonian. The effective dimension 



is the inverse purity of the time averaged state 

1 



d ett (co) 



Tr[w 2 ] ' 



where 



u; = V t = ^e- i (^-s ! )* CfcC r|^)^|. 
k,t 

The Loschmidt echo (see Ref. [30] and the references therein) 
is defined to be the overlap between the initial state and the 
state at time t 



C t = \(i> W t )\ 2 = \(Tp \e- iHt W )\ 2 

= £ e -i(^- Ej )* N 2 W 2. 

k,l 



As 



we have in the non-degenerate case that 

-L=C- t =TPR(iJj )=J2\ck\ 4 . 

k 

Due to the ambiguity in the definition of the eigenbasis, the 
IPR is not well defined in the degenerate case, and 

-jeff = A =^2$E k ,E l \c k \ 2 \c l \ 2 . 
a k,l 



